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We consider weighted random balls in K d distributed according to a random 
Poisson measure with heavy tailed intensity and study the asymptotic behavior of 
the total weight of some configurations in M while we perform a zooming opera- 
tion. The resulting procedure is very rich and several regimes appear in the limit, 
depending on the intensity of the balls, the zooming factor, the tail parameters of 
the radii and the weights. Statistical properties of the limit fields are also evidenced, 
such as isotropy, self- similarity or dependence. One regime is of particular interest 
and yields a-stable stationary isotropic self-similar generalized random fields which 
recovers Takenaka fields, Telecom process or fractional Brownian motion. 
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Introduction 

In this work, we consider the so-called weighted random balls model and investigate 
its convergence when suitably rescaled and normalized. We exhibit three different 
asymptotic regimes driving the macroscopic and microscopic variations of this model, 
namely (i) a stable, translation and rotation invariant, self-similar random field on 
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M. d , (ii) a Poissonian field and (Hi) a stable field with independence. The weighted 
random balls model is constructed in the following way: the centers of the balls 
are distributed according to a Poisson point process, with each center x labelled 
with a random radius r and a random weight m. The field under study is, roughly 
speaking, at each point, the weight density defined as the sum of the weights of 
the balls containing this point. The overlap of the balls yields non-trivial spatial 
correlations when the random radii of the balls are heavy tailed. 

This fairly simple geometric construction has found numerous applications and is 
pertinent in various modeling situations. Similar stochastic models were considered 
by Kaj in [6] when modeling a simplified wireless network that consists of a collec- 
tion of spatially distributed stations equipped with emitters for transmission over 
a common communication channel. Here, the location of a station or of a network 
node is represented by the point x, its range by the radius r and its power by the 
weight m. The weight density measures the total power of emission at a given point 
and in this case, m is supposed to be non-negative. But our model supports more 
generally real- valued weights. 

In p], Bierme and Estrade consider similar models in dimension d = 2 as models 
in imagery (in this case, the weight intensity stands for the gray level of a pixel in 
a black and white picture) and in dimension d = 3 for modeling three dimensional 
porous or heterogeneous media (here, the weight density is seen as a mass density). 
They investigate the microscopic properties of the random balls configurations by 
performing a scaling operation which amounts to zoom in smaller regions of space. 
In [7], Kaj et al. study similar random grain model by shrinking to zero the volume 
of the grains. This amounts to analyse the macroscopic properties of the random 
balls configurations by performing a scaling operation which amounts here to zoom 
out over larger areas. 

Recently, Bierme, Estrade and Kaj introduce in [2] a general framework for 
rescaled random balls model allowing both zoom-in (as in p]) and zoom-out (as in 
[7]). In this zooming procedure, several limit fields arise, which are either of Gaussian 
or of Poisson type according to the respective asymptotic of the zooming rate and 
of the Poisson intensity of the balls. Furthermore, they show that essentially all 
Gaussian, translation and rotation invariant self-similar generalized random fields 
can be obtained as such a limit. 

Note that in the rescaled random balls model of p], [7] and [2], the weights in 
the field under study are fixed equal to m = 1. Models with randomized weights 
have been less intensively studied. In dimension d = 1, Kaj and Taqqu study in [4] 
limiting schemes for weighted random balls model, deriving Gaussian, Poisson and 
stable regimes. This model applies in particular to study the random variation in 
packet networks computer traffic. 

Our main contribution in this paper is to introduce a general study of macroscopic 
and microscopic variations in weighted models in M. d . This generalizes both [2] since 
the balls are randomly weighted and [4] since we consider an arbitrary dimension 
d and more general configurations on the balls. As in [7] and [4], three different 
regimes appear according to the relative behavior of the scaling rate and of the 



Poisson intensity. In particular, when the random weights are heavy tailed, the 
limit generalized random fields are stable, translation and rotation invariant, and 
also self-similar. The paper is organized as follows. The model under study is 
described in Section [U Our main results under different scaling regimes are stated 
and discussed in Section El Finally, Section [3] is devoted to the proof of technical 
lemmas and of the main results. 



1 Model of weighted random balls 

We consider random balls B(x, r) = {y E R d : \\y—x\\ < r} with weight m, the triplet 
(x,r,m) being distributed according to a Poisson random measure N\(dx,dr,dm) 
onl^xf+xR with intensity 

n(dx, dr, dm) = XdxF(dr)G(dm) 

where A is positive, F is a positive measure on R + and G a probability measure on 
R. Here, and in what follows, || • || stands for the usual Euclidean norm on R . 

The point process of the centers of the balls in M. d is the projection of the point 
process in R xR + xR corresponding to the Poisson random measure N\(dx, dr, dm). 
It is easily seen that it is a Poisson point process with intensity Xdx, and hence the 
parameter A is interpreted as the intensity of the balls in R rf . 

We suppose that the measure F driving the distribution of the radius r is abso- 
lutely continuous F(dr) = f(r)dr with 



r a F(dr) < +00 (1) 

and such that for either e = +lore = — 1, 

f(r) ~ r _ e Cpr~ l -P (2) 

where by convention +1 = and _1 = +00. As will be explained later, the 
case e = +1 will be referred as the zoom-in case, whereas the case e = — 1 will be 
referred as the zoom-out case. Condition ([2]) assumes a power behavior of the radius 
density at the origin (zoom-in case e = +1) or at infinity (zoom-out case e = —1). 
Condition (|TJ) is equivalent to the finiteness of the volume of the random balls. Note 
that assumptions (Q]) and (J2]) together imply that for e = +1, we must have (3 < d, 
while for e = —1, we must have > d. 

We suppose that the probability measure G belongs to the normal domain of 
attraction of the a-stable distribution S a (a, b, r) with a £ (1, 2], i.e. if X±, . . . , X n 
are independent and identically distributed (i.i.d.) according to G, rT l l a {X\ + • • • + 
Xn) => S a (a,b,T). We recall the following estimate (see [3]) of the characteristic 
function ipc of G as 9 — > 

<p G (0) = l + i6r- a a \6\ a (l + ibe(9) tan(vra/2) + o(|0| a ), (3) 



where here, and in what follows, e(a) = +1 if a > 0, e(a) = — 1 if a < and 
e(0) = 0. In case a G (1,2), typical choices for G are heavy tailed distributions 
while for a = 2, G is any distribution with finite variance. In this latter case, we 
recover a weighted version of the main results in [2] (set G = 6± to recover exactely 
the setting described in [2]). 

Let M. denote the set of signed measures on M. d with finite total variation |/i|(R ), 
where \fj,\ is the total variation of a measure /i. We recall that equipped with the 
norm of total variation ||ju||x = |jU.|(M d ), M. is a Banach space. We consider the 
random field 



M(n)= mp(B(x,r))N x (dx,dr,dm) (4) 

./R d xR+xR 

indexed by signed measures p G M.. When p, = 6 y , y € M rf , M(5 y ) is the weight 
density at point y as described in the introduction: it is the sum of the algebraic 
weights of the balls containing the point y. 

Note that the stochastic integral in ([4]) is well defined and has finite expected 
value since 



< 



\mp(B(x,r))\n(dx,dr, dm) 

<M+xR 

f \m\G(dm) x \\B(0,l)\\n\(R d ) [ r d F(dr) < +oo. 



where \A\ stands for the Lebesgue measure of a Borel set A. Furthermore, the 
expected value is given by 

E[M(/x)] = A|fl(0,l)| / mG{dm) I r d F(dr) fi(R d ). 



We are interested in the variations of M{fj) at a microscopic or macroscopic level. 
To do so, we swell, resp. shrink, the volume of the balls replacing the radius r of a 
ball by pr and taking the limit p — ► +oo, resp. p — > 0. In this procedure, the law of 
the radius is replaced by F p (dr) = f(r/p)dr/p, the image measure of F(dr) by the 
change of scale r h-> pr. In order to derive non-trivial asymptotics, the intensity A 
of the balls is changed accordingly and we shall write X(p) to underline that from 
now on the intensity depends on the scaling parameter p. In what follows, we are 
thus interested in the following random field: 

M p (n)= I mp(B(x,r))N X ( p ) } p(dx,dr,dm) 

JR d xR+xR 

where Nx( p \ p (dx, dr, dm) is the Poisson measure with intensity X(p)dxF p (dr)G(dm) . 
The limit p — > is interpreted as zoom-out in the random configurations of balls and 
this is relevant when the behavior of / is known at +oo, i.e. e = — 1 in ([2]). In this 
case, we investigate the macroscopic variations of M. On the contrary, p — » +oo is 
interpreted as zoom-in in space and this is relevant when the behavior of / is known 
at 0, i.e. e = +1 in ([2]) and this is the microscopic variations that are investigated. 



Remark 1.1 As observed before, the choice G = 5\ recovers the setting of [2j for 
non- weighted random balls, see (4) therein. If d = 1, a verbatim replacement of 
B(x, r) = (x — r,x + r) by (x, x + r) and the choice /x = | • n(0, t)\ recover the field 
studied in [I] in the "continuous flow reward model", see (18) therein. 

2 Results 

We exhibit normalization terms n{p) such that the normalized centered random 
field n(p)~ l [M p (-) — E[M p (-)]) converges in finite- dimensional distribution (f.d.d.) 
to a limit random field. In what follows, we are interested in f.d.d. convergence on 

subspaces M. of M. and we will denote it by — >. 

It is natural to investigate first the behavior of the random field giving the density 
of the weights at each point which in our notations rewrites (M p (5 y )) yGK d. The 
heuristic is the following. The average numbers of balls containing the point y is 
given by 

E / l{yeB(x,r)}N\(p) tP (dx,dr,dm) = VX{p)p d , 

jR d xK+xM. 

where V = Cd f r d F(dr) is the expected volume of a random ball and c<j stands 
for the volume of the Euclidean unit ball in M. d . Since the weights belong to the 
domain of attraction of an a-stable distribution, it is natural to introduce the scaling 
uq(p) = X(p) 1 ' a p d ' a . Convergence of the normalized and centered random variable 
M p (5y) to an a-stable distribution is obtained if we suppose that \{p)p d — > +oo 
when p — » -e . Heuristically, the dependence between M p (5 yi ) and M p (5 y2 ) is given 
by the weights of the balls containing both points y\ and ?/2- In the zoom-in case 
(e = — l,p — ► +oo), the balls are very large yielding total dependence at the limit 
and we have: 

n { P y l (X p (5 y ) - E[X p {S y )]) f -^W a , y G R d (5) 

where W a (y) = W a is a constant random field distributed according to S a (aV 1 ' a , b, 0). 
In the zoom-out case (e = — 1, p — > 0), the balls are very small yielding independence 
at the limit and we have: 

no{p)-\M p [8 v ) - E[M P {S V )]) ^ W a (5 y ), y e R d , (6) 

where W a (S y ),y € M. d , are i.i.d. S a (aV 1 ' a , b, 0) distributed. Similar results as in 
(EJ) and in ([6]) hold true for f.d.d. convergence on the space of measures with finite 
support. Since these results are not surprising, their proofs are omitted and in what 
follows we investigate convergence for more general measures. 

2.1 Preliminaries on measured spaces 

We introduce a subspace M a ,p C M. on which we will show the convergence of the 
rescaled generalized random field M p (p). 



Definition 2.1 For 1 < a < 2 and (3 > 0, let M a ,p be the subset of measures 
fi G M satisfying for some finite constant C and some < p < (3 < q: 

7(r) := / \n(B(x,r))\ a dx<C(r p Ar q ) (7) 

JR d 
where for reals a,b: a f\b = min(a, b). 

Here, and in what follows, C is a finite constant that may change at each oc- 
currence. Some elementary properties of the spaces A4 a ,/3 are given in the following 
proposition. 

Proposition 2.2 

i) M a ,p is a linear subspace of M. on which 



\/fi G M a 



p , / \n(B(x, r^fr-P^dxdr < +oo. 

jR d xR+ 



ii) M. a fi is closed under translations, rotations and dilatations, i.e. when /i € 
•M a ,0, Tsfj,, Qfi and \x a are also in M a ,p where for any Borelian set A and for 

sel d ,0e 0(R d ), ael+ 



T s fi(A) = fi(A - s), Ofi(A) = f i(e- 1 A),^ a (A) = ii{a- l A). 

Hi) When a < a' , we have A4 a p C M a > p. 

iv) When j3 > d, the space M a ,f3 * s included in the subspace of diffuse measures 
(i.e. such that /i({x}) = for any x € M. d ). 

v) When j3 < d, the space M a ,p is included in the subspace of centered measures 
(i.e. such that fi(R d ) =0). 

Observe that Dirac measures 5 y , y G M. d , are not in M. a ^. However, explicit ex- 
amples of measure in M a ,p are given in the following proposition. Absolutely con- 
tinuous measures (with respect to the Lebesgue measure) fi(dx) = (j)(x)dx with 
integrable density <f> € L 1 (M d ) nL a (M. d ) will play an important role. In this case, we 
shall (abusively) note \i G L l (R d ) n L a (R d ). 

Proposition 2.3 

i) If d< (3 < ad, any measure fi G L 1 ^) DL a (M. d ) belongs to M a ,p- 

ii) If d — 1 < j3 < d, any centered measure fi(dx) = (p(x)dx G L 1 (IR c( ) such that 
L d ||y|||</>(y)|dj/ < +oo belongs to M a ,p, as well as any centered measure with 
finite support. 

Note that in particular, when d < (3 < ad (resp. d — 1 < f3 < d), M. a ,p contains 
the space S of measures with density in the Schwartz class (resp. Sq the space of 
centered measures with density in the Schwartz class). Note also that when a = 2, 
the conditions supposed in [2] on the measure fj, (expressed in terms of Riesz energy) 
imply that \i G Mifi- By analogy with the case a = 2, we suspect the space M. a ,p 
to be reduced to {0} whenever (3 < d — lor/3> ad, but we have no formal proof of 
these facts. However, we refer to Theorem 12.191 for a positive result when (3 > ad. 



2.2 Limit theorems for the rescaled weighted random 
balls model 

We now come to the main results of this paper, viz. limit theorems for the rescaled 
generalized random fields M p and for configurations p € M. a ,p on the balls. As in [7] 
and [4] (for e = — 1), several regimes appear according to the density of large/small 
balls in the limit. More precisely, using ((2)): 

Zoom-out case (e = —1, i.e. (5 > d and p — ► 0). The mean number of balls 
with radius larger than one that cover the origin is given by 

/ l M<r l r>1 X(p)dxF p (dr) = c d X(p) [ ™r d F p (dr) ~^ ^A(p)/. 

jR d xK+ J\ p — a 

Consequently, we distinguish the following three scaling regimes: 

• large-balls scaling: X(p)p@ —* +oo, 

• intermediate scaling: X(p)p^ — > a G (0,+oo), 

• small-balls scaling: X(p)p 13 — > 0. 

Zoom-in case (e = +1, i.e. (3 < d and p — > +oo). The mean number of balls 
with radius less than one that cover the origin is given by 

/ l M<r l r<1 X(p)dxF p (dr) = c d X(p) [ r d F p (dr) ~ p ^ +00 ^|A(p)/. 

JR d xiR+ Jo a — p 

In this case, the three scaling regimes are: 

• small-balls scaling: X(p)p 13 — > +oo, 

• intermediate scaling: X(p)p@ — > a G (0,+oo), 

• large-balls scaling: X(p)p@ — > 0. 

In what follows, we study precisely the limiting shape of the random balls by 
investigating the fluctuations of M(p) around its mean. Three different limit fields 
are exhibited according to the scaling performed. 

2.2.1 Stable regime with dependence 

In this section, we investigate the behavior of M under the scaling p f3 X(p) — > +oo. In 
this case the limiting field is given by an a-stable integral. We recall that the stable 
stochastic integral of / with respect to an a-stable random measure with control 
measure m is well defined whenever / G L a (dm) and in this case, this stochastic 
integral follows an a-stable distribution. We refer to [10] for a complete account on 
stable measures and integrals. The asymptotic of the rescaled generalized fields M p 
is given by the following result: 



Theorem 2.4 Suppose p^A(p) — > +00 when p — » e . Lei ni(p) = \(p) 1 ' a p^' a . 

We have 

M P H-E[M, ( .)j^ _, 

where Z a (p) = / K d x jg+ p(B(x,r))M a (dr,dx) is a stable integral with respect to the 
a-stable measure M a with control measure a a Cpr~ ~"drdx and constant skewness 
function b given in the domain of attraction of G. 

Note that Z a (p) makes sense as soon as f^_ x ^+ \p(B(x,r))\ a r~ 1 ~^drdx < +00 (see 
Prop. 12.21 -i)). However, we need the stronger assumption p G M. a ,p m order to 
derive (|8|). Roughly speaking, the control (j7]) of p G M. a ,/3 allows to replace F by 
its tails behavior given in (J2j) in asymptotic estimate. 

Due to the invariance by translation and rotation of the Lebesgue measure, the 
self- similarity of stable integral and the (global) invariance by rotation of the balls 
and because of Proposition I2.21 ii). we derive the following properties for the limit 
field Z a of Theorem E3J 

Proposition 2.5 

i) The field Z a is stationary on M a ,p, that is: 

Vp € M a , p , Vs G R d , Z a (r s p) f = Z a {p). 

ii) The field Z a is isotropic on M a ,p> that is: 

V/i G M a ,p, V6 G 0(R d ), Z a (Qp) f = Z a {p). 

Hi) The field Z a is self-similar on M. a fi with index [d — j3)/a, that is: 

\/p G M a ,p,\/a > 0, Z a (p a ) f = a^ d -^/ a Z a {p). 

Remark 2.6 The covariation gives an insight into the structure of the spatial de- 
pendence of the stable generalized field. It is a generalization of the usual notion 
of covariance to the stable framework. Here, for p\,p<i G M a ,p, the covariation of 
Z a (pi) on Z a (p 2 ) is given by 

[Z a (pi),Z a (p 2 )}a 

= ° a Cp f p 1 (B(x,r))e(p 2 (B(x,r)))\p 2 (B(x,r))\ a - 1 r-f 3 - 1 drdx. 

JR d xR+ 

Note that the integral above is well defined by Holder's inequality since p\ and p 2 
belong to M. a ,/3- We refer to [10] for a definition and properties of the covariation. 
Note that unlike the Gaussian case, the covariation structure is not sufficient to 
characterize the distribution of the generalized random field. However, since even 
if pi and p 2 have disjoint supports, [Z a (pi), Z a (p 2 )] a ^ 0, Z a {p{) and Z a (p 2 ) are 
not independent and the random field Z a is stable with dependence. 

S 



Remark 2.7 Note that when d — 1 < (3 < d, p z = 5 Z — do for z € R d belongs 
to M a , (3- For such a measure, when moreover 6 = (i.e. when G in our model is 
symmetric), our limiting field rewrites 

Z a {Hz) = / 1 B(z,r)AB(0,r) M a(dx,dr) 

JR d xR+ 

where AAB = (A \ B) U (B \ A). In this case, we recover the so-called (a,H)- 
Takenaka field with H = (d — (3) /a. It is self-similar with index if, with stationary 
increments and almost surely with continuous sample paths, see [2j p. 25] or pjH 

Sect. 8.4]. 

Remark 2.8 When d = 1, f3 & (1, a) and /it = \ ■ fl(0, t)\, the field Z a (fj, t ) coincides 
with the Telecom process obtained in the fast connection rate for the "continuous 
flow reward model" in [U Th. 2], see also Remark ll.ll above. Moreover for a = 2, 
-^(a**) is a fractional Brownian motion of Hurst index H = (3 — j3)/2 E (1/2,1) 
(note that, for a > 0, ^ a t(A) = a[i t (a~ 1 A)). 

Remark 2.9 When a = 2, Theorem 12.41 exhibits a Gaussian limit field and gener- 
alizes Theorem 2.1 in [2] with random weights. Indeed, in this case, we have (up to 
some multiplicative constant) Z2 = Wp. 

Remark 2.10 A natural complementary result to be investigated is the tightness of 
M p after normalization and centering which would allow to turn f.d.d. convergences 
into weak functional convergences. In dimension d = 1, only partial tightness results 
are available for the processes studied in [4], [9] (see Section 4 on "continuous flow 
reward model" in [4j and the remarks of Th. 1, Th. 2 and Th. 3 in |9j). In the 
case of generalized random fields, tightness issue is more difficult to tackle due to 
the lack of tractable tightness criterion. 

2.2.2 Poissonian regime 

In this section, we investigate the behavior of M under the scaling p^\(p) — ► a 6 
(0,oo). In this case, the limiting field is given by a compensated Poisson integral 
and we refer to [8] for a general description of Poisson integral. We have: 

Theorem 2.11 Suppose X(p)p^ — > a d ~@ when p — » 0~ e for some a > 0. We have 

M p {p) - E[M p (jj)] — -> J(p a ), P -» (T e 
where /i a is the dilatation of p and J is the compensated Poisson integral 

Jil 1 ) = / mp(B(x,r))Np(dx,dr,dm) (9) 



with respect to the compensated Poisson random measure Np with intensity given by 
Cf3r~P~ 1 dxdrG(dm) . 

9 



Note that the Poisson integral in ((9)) above is well defined since 

(\mp(B(x,r))\ A (mp(B(x,r))) 2 ) r' f3 ' l dxdrG(dm) < +00 (10) 



see Section l3~4l As the stable field Z a , the Poisson field J enjoys similar properties. 
However, note that in contrast to Z a , J is not self-similar but (and similarly to [2], 
see also [5]) J satisfies an aggregate similarity property. 

Proposition 2.12 The field J is stationary and isotropic on M a ,p- Moreover, J 
is aggregate similar, viz. V// £ M. a ^,\/m > 1, 



j(vaj f = xyo.) (11) 



fdd- 

where J 1 , 1 < i < m, are independent copies of J and a m = m • ^ ~"' . 



The proof of this proposition follows from straightforward computation and will be 
omitted. A comparison of the limiting procedures in Theorem 12.41 where \(p)fr — > 
+00 and in Theorem 12.111 where X(p)p^ — * a d ~@ suggests that when a d ~P — ► +00, we 
can recover Z a from J . This is true and precisely stated in the following proposition: 



Proposition 2.13 When a d @ — > +cxo, we have , d } )3) i a J(p a ) — ■* Z a (p). 



Remark 2.14 As in Remark 12.81 when d = 1 and \xt = \ ■ fl(0, t)\, the field J(p-t) 
coincides with the intermediate Telecom process obtained in the intermediate con- 
nection rate for the "continuous flow reward model" in [4j Th. 1], see also Remark 
1.11 above. 



Remark 2.15 When a = 2, Theorem 12.111 generalizes Theorem 2.5 in [2] with 
random weights. The field J recovers Jp in [2] when the random weights in our 
model are constant. Otherwise the law of J depends on the law G of the weight. 

2.2.3 Stable regime with independence for small radius 

In this section, we investigate the behavior of M under the scaling p^X(p) — ► 0, 
but we restrict to the case d < f3 < ad, i.e. e = — 1 and p — > 0, see Section [2.2.41 
for (3 > ad. The case m = 1 is considered in Theorem 2 Hi) of [7] and we extend 
here the results and proofs to the case when the weights are random and belong to 
the normal domain of attraction of a stable distribution. In comparison to the case 
(3 < d, the tails of the law of the radius are lighter and thus the radius considered 
are small. We show that the asymptotic behavior is given again by a stable field but 
with index 7 = /3/d and defined on M. d . Moreover in contrast to the stable field Z a 
of Section [2.2.11 this new field exhibits independence. 



10 



Theorem 2.16 Let d < (3 < ad and suppose that A(p) — > +00 and \{p)p^ — > as 
p — > 0. Then with ri2(p) ■= \(p) '@p and 7 = (3/d S (l,a), we have 

M p (-)-E[M p (-)} LHR d )nL<*(R d ) 

Mp) 

where, for p{dx) = <p(x)dx, Z^(p) = f Rd 4>(x)M^(dx) is a stable integral with respect 
to the "{-stable measure M 7 with control measure Oydx for 

4-^f i^*/|mrG( dm ) 

a JR + r ' JR 

and with constant skewness function equals to 

f R e(m)\mpG(dm) 
b i- / R |m|7G(dm) • (U) 

Note that the integrals above are well denned when d < [3 < ad (see Lemma 13.11 
below). The limiting field Z 7 enjoys similar properties as Z a and J: 

Proposition 2.17 The field Z^ is stationary, isotropic and self-similar with index 

(d-/3)h. 

Remark 2.18 As in Remarks 12.81 and 12.141 when d = 1 and <\>t = l(o,t)j the field 
Z-y(4>t) coincides with the process obtained in the slow connection rate for the "con- 
tinuous flow reward model" in [H Th. 3], see also Remark 11.11 above. In this 
particular case, Z^(4>t) is a 7-stable Levy process. 

2.2.4 Stable regime with independence for very small radius 

When the tails of the radii are lighter than that in Section 12.2.31 i.e. (3 > ad, 
the same stable regime with independence as in Section 12.2.31 appears but under 
a different normalization n-${p) := \{p) l ' a p d and a different stability index a. As 
previously, since (3 > ad, we have e = — 1 and the limits are taken when p — ► 0, i.e. 
the limiting scheme is a zooming-out procedure. 

Theorem 2.19 Let (3 > ad and suppose that X(p) — > +00 as p — > 0. Let n^(p) := 
\{p) l / a p d , then 

Mp{-)-E[M p {-)} L\R d )nL<*(R d ) 

where, for p(dx) = (j)(x)dx, Z a (p) = J^ d cj)(x)M a (dx) is a stable integral with respect 

to the a-stable measure M a with control measure o~ a dx with a a = acd ( L + r ad F{dr)) 
and constant skewness equal to b. 
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Remark 2.20 It is worth noting that in both Theorem 12.161 and Theorem 12. 191 the 
stable regime is driven by the parameter 7 = (f3/d) A a, since the normalization is 
A( / o) 1 ' 7 p d and the stability index is 7. 

Actually, only the asymptotics of the law with the heavier tails contribute to 
the limit while the law with the lighter tails appears only (but globally) as a mere 
parameter in the limit. In particular, observe that Theorem 12.191 applies for any 
distribution F such that L + r ad F{dr) < +00. 

Remark 2.21 When d = 1 and pt = \ ■ n(0, t)\, we recover ii) in Theorem 4 of [I]. 

3 Proof of the results 

In what follows, note that the linearity of the random functionals M p and of the 
stochastic integrals in W a , W a , Z a , J and Z 7 , together with the Cramer- Wold device 
imply that the convergence of the finite-dimensional distributions of the centered and 
renormalized version of M p is equivalent to the convergence of the one-dimensional 
distributions. To do so, we will explicitly compute the limits of the characteris- 
tic functions, denoting tpx for the characteristic function of a random variable X. 
Observe that the characteristic function of n(p)~ 1 (M p (p) — K[M p (fi)]) rewrites: 

¥>n(p)-i(M p ( M )-E[M p ( M )])(#) 

= exp ( / ^ {n{p)~ l 0mp{B(x,r))) \{p)dxF p {dr)G(dm) 



where \P(u) = e lu — 1 — iu, see [8]. Integrating first with respect to the probability 
G(dm), we have 

Vn(p)-i(M p (n)-E[M p (n)])(d) 

= exp( [ X(p)^ G (n(p)~ x e^B(x,r)))dxF p (dr)) (13) 

\jR d xR+ J 

where ^g(u) = J R x ^(mu)G(dm). We also recall that the characteristic function of 
the stable distribution S a (a,b,r) is given by exp(—a a \x\ a (l — ibe(8)tan(Tra/2)) + 
ird). 

3.1 Preliminary lemmas 

In this section, we collect some useful lemmas that will be needed in the proof 
of our limit theorems 12.41 12.111 and 12.161 We recall the following estimate for the 
characteristic function of distribution in the domain of attraction of a stable law: 

Lemma 3.1 Suppose X is in the domain of attraction of an a-stable law S a (cr, 6, 0) 
for some a > 1. Then 

cp x (0) - 1 - i9E[X] ~ -o- a \9\ a (l - ie{9) tan(vra/2)6). 
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Furthermore, there is some C > such that for any flel, 

\<p x (O)-l-i0E>[X\\ <C\d\ a . 

The following lemma is a reformulation from lemma 2.4 in [2]. It shows that in 
the scaling limit p — > _e , the behavior of F p is given by the power tail of F. This 
is crucial in several estimates. 

Lemma 3.2 Let F be as in (J2J) and e = ±1. Assume that g is a continuous function 
on R + such that for some < p < (3 < q, there exists some C > such that 

\g{r)\ <C(r p Ar«). (14) 

Assume furthermore that (g p ) p >o is a family of continuous functions such that 

lim \g(r) — g p (r)\ = and \g(r) — g p (r)\ < C(r p A r q ). (15) 

Then 

g p (r)F p (dr) ~ C^p^ / g{r)r~ l ~ l3 dr when p — > CT 6 . 

In the proof of Theorem 12.41 and of Theorem 12. Ill below, this lemma will be used in 
the particular case where g p = g and g satisfies condition (JT4J) . Roughly speaking, 
the proof of Lemma E21 consists in taking the limit in the integral. This is authorized 
by the dominated convergence theorem under (fi"4"l) and (JT5J) . We refer to [2] for more 
details. 

3.2 Proofs of Propositions E21 and [231 

Proof of Proposition 12.21 Proof of i). If (JT]) holds true for p\ with p\ < f3 < q\ 

and for p,2 with P2 < f3 < q2, then ([7]) holds true for p>\ and p,2 with p = p\ V P2 < ft 
and q = qi A q2 > fl (possibly with a different constant C). For all a\,a2 € R: 

/ \(aip 1 +a 2 p2)(B(x,r))\ a dx = \\(aipi + a 2 p2){B(x,r))\\% 

< (|oi|||Aii(B(z,r))|| Q + \a2\\\p2(B{x,r))\\ a ) a 

< ((\ ai \ a C(r p A r 9 )) 1 /" + (\a 2 \ a C(r p A r^f^Y 
= C(\a l \ + \a2\) a (r p Ar q ). 

This is ([7} for a\p\ + 0,2^2- 

Proof of ii). Since (r a fi)(B(x,r)) = p(B(x-s,r)), (6p)(B(x,r)) = ^(e^x^)), 
p, a (B(x,r)) = /i(S(a _1 x,a _1 r)), the closeness of M a ,p by translations t s , by rota- 
tions G and by dilatations x \— ► ax follow straightforwardly from the invariance of 
the Lebesgue measure by translations, by rotation, and by an immediate change of 
variable in ([7]). 
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Proof of iii). Since |/j|(M d ) < +00, for [i G M a ,f5 and a < a' , we have 
\fi(B(x,rW'dx [ \n(B(x, r))\ a '- a \n(B(x,r))\ a dx 

\n{B(x,r))\ a dx 



< \fi\(R d ) a - c 

< C(r p A r q ) 



which proves \i G Ai a 



'fi- 



Proof of iv). We prove that \x G Ai a ,0 is diffuse when j3 > d. Indeed, suppose 
that /j, has an atom a, then for small enough r, 7(7") > \fi(a)/2\ a c,ir d , where we 
recall that 7(7") is defined in ([7]). To see this, let e > be such that ||/i|(5(a, e)) — 
|/x(a)|| < |/x(a)|/2. Then, for every r < e/2 and x G B(a,r), \fi(B(x,r))\ > \/j,(a)\/2. 
Integrating on x G B(a,r), we get 7(7*) > (|^(a)|/2) a c,ir d . This is in contradiction 
with ([7|) which rewrites 7(f) < Cr 9 for q > (3 > d when r is small. 

Proof of v). We prove that \x G M a ,p is centered when (3 < d. We will show that 

7 (r) > |/x(M d )/3| a Cd r d (16) 

when r is large enough. This is sufficient to prove v) since (jj|) rewrites 7(f) < Cr p 
for p < P < d when r > 1 which is in contradiction with (fi~6l) when /i(M d ) 7^ 0. 

The bound fEJ is obvious if /i(M d ) = and if fi(R d ) ^ 0, let M be such that 
\H\(B(0,M) C ) < \n(R d )\/3. Then, for r > M and any x G 5(0, r - M), B(0,M) C 
B(x,r) and |^(S(x,r))| > |/i(M d )| - \n\{B(x,r) c ) > 2|/x(M d )|/3. Integrating on 
x G 5(0,r - M), we obtain 7(7-) > (2|^(R d )|/3) a c d (r - M) d . This implies (fT6jl . 



Proof of Proposition 12.31 Proof of i). First, when d < j3 < ad and n{dx) 

(p(x)dx G L 1 ^) n L a (R d ), we have: 



/ 



\fi(B(x,r))\ a dx 



-{\\x—y\\<r 



}4>(y)dy 



dx 



< 



l {\\x-y\\<r}\<f>(y)\ a dV 



Q-l 

l{\\x-y\\<r}dy I d.r 



{c d r d r~ l 

(c d r d ) a 



[ {\\x-y\\<r}\Hy)\ a dydx 



l {\\x-v\\ 



\\<r}dx) \(p(y)\ a dy 



\Hy)\ a dy 
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where we applied Holder's inequality with a > 1. Next, in the same way, 



/ \iM(B(x,r))\ a dx 

JR d 



< 



dx 

a-l 



i{\\x-y\\<r}<i>{y)dy 
U\\x-y\\<T\\<l>{v)\dv) ( / \4>{y)\dy) dx 

i ' J \Jm. d J 

c d r d ( [ \4>(y)\dy] . (17) 



As a consequence, condition ([7J holds with p = d<(3<q = ad, and /i 6 M a ,p- 

Proof of ii). Suppose d — 1 < (3 < d and /j,(dx) = 4>(x)dx e L^R^) is centered. 
Using n(M. d ) = 0, we have: 

/ \n(B(x,r))\ a dx (18) 

JR d 

= (liU-vHr} ~ 1 {\\x\\<r}) < t>{y)dy dx 

Jm. d JR d 

- / (/ l 1 {iix-i/ii<r} - 1 {||x|]<r}ri<? i (y)l^) ( / \4>(y)\dy) dx. 

JR d \JR d J \JR d J 



Let 0{z) = |-B(0, l)AB(z, 1)| denotes the volume of the symmetric difference of the 
balls with unit radius centered at and at z £ R d . We have, 



da(V 



/ iMWx-y^r} ~ l{\\x\\<r}\dx = V 

Jx£R d Kr 

The function 6 is continuous, upper bounded by Cd and 9{z) = 0(||z||) as z — ► 0. As 
a consequence, the global estimate \0(z)\ < C\\z\\ holds true for some C > 0. This 
entails 

| A *(S(x J r))| a d a: = ( f \cf>(y)\dy] i r d 8 (?-) <j>(y)\dy (19) 

\jR d J JR d ^rJ 



a-l 



< / \<Ky)\dy) / C||y|||^(y)|di/r d - 1 

\J]R d / JR d 

< CV* -1 . 

As a consequence, condition holds true with p = d — l<(3<q = d because (fTTl) 
still holds true, and finally \x £ M. a ,0- 

Alternatively, if \x has a finite support {oi,. . . ,%>}, let (5 > such that for 
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1 < i < p, B(a,i,5) n Supp (/*) = {ai}. For r < 5/2, 

7 (r) = I \p(B(x,r))\ a dx 

= JZJ HB(x,r))\ a dx 

= c d ^T\p( ai )\ a r d = 0(r d ). (20) 



i=l 



Next, let M be such that /i(5(0,M) c ) = and note that p(B(x,r)) = when 
B(x,r) n 5(0, M) = or when 5(0, Af) C S(x,r) since /x(R d ) = 0. We derive 
/i(5(x,r)) = when ||x|| < r — M or when ||x|| > M + r. Since ^ is a finite measure, 
we have 



7 (r) = / H5(£,r))| Q dx 

Jr-Af<||x||<r+Af 

< c d ((r + M) a! -(r-Af) d )(|^|(IR d )) Q 
= 0(r d - 1 ), r^+oo. 

Together with (|20l) . this yields condition ([7]) with p = d — 1 < j3 and q = d > f3. D 



Remark 3.3 (On the bound for large radii) Note that in order to derive the 
bound 7(7*) < r p for p < j3 when r is large, the existence of a density for \x is 
not required. We can instead suppose that \x satisfy some tail condition: for some 
fj > d/a 

|/i|(5(0, R) c ) = 0{R~ f <) as R -» +00. 

3.3 Proof of Theorem [231 

The characteristic function of the stable integral Z a (/j,) is given by 

<pz a M(e) (21) 

= exp \-C p a a f \0p(B(x, r))|°(l - ie(9fi(B(x, r))) tan(7ra/2)&)r- 1-/3 dr(ia; J . 

Since the characteristic function of the Poisson integral n\(p)~ l (M p {p) — K[M p (fi)]) 
is given by (PT3]) . comparing (|2T1) and (Tl3l) . it is sufficient to show that 



lim / \(p)^ G {ni{p)- 1 eix{B{x,r)))dxF p (dr) (22) 

-Cpcr a I \0n{B(x, r))\ a {l - ie(6p(B(x, r))) ttaifrafflfyr^-Pdrdx. 

Jm. d xR+ 
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Since n\(p) = (X(p)p^) 1 ' a — > +00, Lemma l3~Tl applies and yields 
\{p)^ G { ni {p)- l ep{B(x,r))) 



f~\_l 



a a p~^\6\ a \pL(B{x,r))\ a (l- ie{6p(B(x, r))) tan(vra/2)6). 



Since | n , s /j,(B(x, r))| < - ? -t |/i|(R ), this equivalence relation is uniform both in 2 
and r and can be integrated. This yields 

A(p)* G {nxipT^niBix, r))) dxF p {dr) (23) 

( xR+ 

~ -a a p- p \9\ a I \p(B(x,r))\ a (l-ie(6p(B(x,r)))tan(ira/2)b)dxF p (dr). 

jR d xR+ 
Finally Lemma 0T2] applies with 

g(r)= / \p{B(x,r))\ a (l-ie{6p{B(x,r)))tan(Tra/2)b)dx, 

jR d 

note that ([7]) implies that 5 satisfies condition dH|. Consequently 

|/i(5(x,r))| Q (l-ie(6» / u( J B(x,r)))t a n(7ra/2)6)dxF p (dr) 



/ 



xR+ 



~ C^/ / \p J (B(x,r))\ a (l-ie{9p J {B(x,r)))t&n(TTa/2)b)r- /3 - 1 dxdr.(24) 

JR d xR+ 

Finally (j23j) and (|24l) together imply (j22j) . and as explained at the beginning of 
Section El Theorem E3J □ 

3.4 Proof of Condition (HDD 



We prove that Condition ([10]) for the existence of J is satisfied. Note that this 
condition splits into: 

/ (mp(B(x,r))) 2 r'^~ 1 dxdrG(dm) < +00 (25) 

V|m|t(B(x,r))|<l 

and 

\mfi(B(x,r))\r~ l3 ~ 1 dxdrG(dm) < +00. (26) 



J\rr 



\m/i(B (x ,r))\>l 

We shall use the following Lemma for the truncated moments of a distribution in 
the normal domain attraction of a stable law: 

Lemma 3.4 Let G be in the normal domain attraction of an a-stable law for a > 1. 
There are C\,Ci G (0, +00) such that for all x > 0: 

f \m\G(dm) < Cix 1 ^ and f m 2 G(dm) < C 2 x 2 ~ a . 

J\m\>x J —x 
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Proof of Lemma 13.41 Prom [31 XVII. 5], we have f^ x m 2 G(dm) ~ Cx 2 ~ a when 
x — > +oo (note that since G is in the normal domain of attraction, there is no slowly 
varying function in this estimate). But since moreover for x € [0, 1] 

o 
X fX 



m G(dm) = / G{m : u < m < x )du < x < x 

■ Jo 

and the second part is proved. 

Next, since lim^^o i| m i >a . \m\G(dm) = J* R \m\G(dm) < +oo while x 1_a — > +oo, 
x — > 0, the first part comes from [H Eq. (5.21)]: 

2 - a 1 /"^ 2 2 - a x _ a 



□ 



|m|G(iim) ~ / m G{dm) ~ x a , x — > +oo. 



Now, we prove (I25|) and (|26|) . First for (|25l) . we have: 

./|m/j(' J B(2,r, > )|<l 



/• / f +l/\n(B(x,r))\ \ 

< / m 2 G{dm)\ {i(B(x,r)) 2 r-P- l dxdr 

JR d xR+ \J-1/ 



i/HB(x,r))\ 



< C 2 I \n(B(x,r))\ a ' z ^{B(x,r)Yr~^- i dxdr 

1 xR+ 



\ a - 2 li(B(x,r)V--'- L > 

< C 2 I \^{B{x,r))\ a r-^- 1 dxdr 

JR d xR+ 

which is finite when /x S M. a $ (see Prop. 12.21 -i)). Next for (|26j) . we have: 

\mii(B(x,r))\r~"~' dxdrG(dm) 



\mfj,(B(x,r))\>l 

< / \m\G(dm) \n(B(x,r))\r~ ,3 ~ 1 dxdr 

JR d xR+ \J\m\>l/\fi(B(x,r))\ J 

< d [ \ l x{B{x,r))\ a - 1 \n{B{x,r))\r- p ' 1 dxdr 

JR d xR+ 

< d [ \ l x{B(x,r))\ a r~ li - 1 dxdr 

jR d xR+ 

which, again, is finite when fi € M. a ^. □ 

3.5 Proof of Theorem \TU\ 

As in the proof of Theorem [23 it is enough to consider convergence of one-dimensional 
marginals. The characteristic function of the Poisson integral M p (/j,) — E[M p (^)] is 
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given by (fi~3l) and that of the generalized random field </(//) is given by 

(p J(l t a) (0) = expf f *(6'm^(S(a- 1 a; ! o- 1 r)))C / 3r- 1 - /3 drdxG(dm) 

V/R d xR+xM 

= expf/ ^ G (^( J B(x,r)))C /3 a a! - /3 r- 1 - /3 drdx) . 

\JM d xR+ / 

Prom Lemma EH \^ G {9^{B{x,r)))\ < C\6\ a \n(B(x,r))\ a for some C > 0, so that 
condition JT4J for g(r) = L, d tyQ(iJ,(B(x, r)))dx is given again by ([7j when \i € M a ,p- 
Thus, Lemma [3721 applies and together with lim p ^ - e \{p)p 13 = a d_,s entail 



lim / V G (9ii(B(x,r)))dx\(p)F p (dr) 

C p a d -P [ ^GiOfiiBix^Wr-^drdx. 



Since one-dimensional convergence is enough, this achieves the proof of Theorem 
I2TT1 □ 

3.6 Proof of Proposition [27131 

We consider the subsequence a m = rn 1 '^^' . From the aggregate-similarity of the 
field J (see ([PI]) in Proposition 12. 12(1 . we have: 






m 
for independent copies J J , 1 < % < m, of J. But 

= expfm/ 1 $>G{m~ 1 / a 8iJ,(B(x,r)))Cf3r~ 1 ~ l3 drdx ) , 

and from Lemma l3.ll 

y G (m- 1/a 6ti(B(x,r))) ~ a a |#H^(£(x,r))| a (l - ie(0/x(B(a;,r))taii(7ra/2)6). 

The relation above is uniform both in x and r and it is thus integrable with respect 
to drdx. This yields 

= exp ( Cfsa a \6\ a f \p(B(x,r))\ a (l - ie(6 n(B (x , r)) tan ( — ) b) r'^drdx 
V Jm d xR+ ^ V 2 / / 

A standard argument completes the proof of convergence in distribution along an 
arbitrary sequences. □ 
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(27) 



3.7 Proof of Theorem [27161 

We follow the argument in the proof of Theorem 2 in [7j. Recall that here d < f3 < ad 
so that e = — 1 and the limits are taken when p — > 0. Again, by linearity, using the 
Cramer- Wold device, it is enough to deal with one-dimensional marginals. From 
(fj~3j) with a change of variable, the characteristic function rewrites 

= exp ( f * G (en 2 (p)- 1 p(B(x,n 2 (p) 1 / d r))) \(p)dxF {) . l/d {dr\ 

Let p(dz) = cf)(z)dz with e L 1 ^) n L a (R d ), then, from Lemma 4 in [7], as 

«2(p) -> 0, 

n 2 ( / 9)-V( J B(^,™2(p) 1/(i r)) -► 0(x)c d r d 
dx almost everywhere and 

x ,_> 0*(x) = sup (c7 1 tr < V|0B(x,'y)) N ) G L a (M d ). 
■u>o ^ ' 

As a consequence, 

*G ^n 2 (p)-V(S(^,^(p) 1/d r))) A(p)(ixF pn2(p) - 1/d (dr) 

- C7^A(p)p /3 n 2 (p)-^ / \p G (^(x)c d r d ) r^~ x drdx. 

To see this, apply Lemma f3T2l to 

5 (r) = / * G (0<Xx)c d r d ) dx 

and to 

<fc(r)= / ^ G K 2 (p)-V(S(x,n 2 (p) 1 / d r)))dx. 

Since |^g(^)| < C(M A H"), we have 

|«/(r)|<C f min(ca[5|W L ir d ,cg|en|0||2ar ad 

so that condition (fi"4"l) is satisfied with p = d and g = ad. Furthermore, since ^q is 
a if-Lipschitzian function for some finite K , we get 



(2* 



s(r)- 5p (r)|<tfc d r d |0| / c, V^n,^- 1 ^^^^) 1 /^)) - 0(x) 



dx. 



The integrand |c^7 r _d n 2 (p) _1 / u(B(x,n 2 (/9) 1/c( r)) — ^(x)| converges to zero dx al- 
most everywhere. Since its L Q -norm is bounded by ||</>*||L a + ||<^||z, a ; it is uniformly 
integrable and as a consequence, 



lim 



c d 1 r- d n 2 (p)- 1 p(B(x,n 2 (p) 1 / d r)) - <f>(x) 



dx = 0. 
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On the other hand, since for some C > 0, |^g:(v)| < C|u| a , we obtain 

\g(r)-g p (r)\<C(U*\\ La + U\\ La )r ad . 
Hence, g p satisfy condition (fi~5j) with p = d and q = ad. This proves (J28 



Prom the definition of ^(p), ^{p)p^n2(p)~^' d = 1. Furthermore, by splitting 
the integration over R d into {x G M d : 9(j)(x) > 0} and {x G M d : 9<j){x) < 0} and 
performing a change of variable, we have 

I ^ G (e<j){x)c d r d ) r~ p ~ x drdx = D [ (9(j){x)y + dx + D [ (0<f>(x))ldx, 

jR d xR+ ^ ' JR d JR d 

where D is the complex conjugate of D = d^ x cj / R+ \l/G(r)r~ 7-1 dr. We deduce 
t Pn 2 ( P )-HM p (^E[MM]){0) = exp (-o-]|6>| 7 [1 + ib,pe(9) tan (— ^ 



where 

and 

tan(7r7/2))/ R+ r-i-7(l-cos(r))dr / R |m|7G(dm) ^ / Rd Mz)|Tete 

(29) 
But since for 7 G (1, 2), 



f 

.70 



-dx =- \up- — (cos(7T7/2) — ze(tt) sin(7T7/2)) 



x 1+ t ' ' ' (l- 7 )(2-7) 

see Lemma 2 in (3J XVII. 4] (with p = l,q = therein), the first ratio on the right- 
hand side (|29"1) is —1 and we have b^ = bj^-j — , ,, ., 7d — where 6 7 is given in (PT2J) . 
This achieves the proof of Theorem 12. 161 □ 

3.8 Proof of Theorem [2391 



The argument uses the same tools as in the proof of Theorem 12. 161 and we only give 
here the main lines. From (fi~3j) and a change of variable, the characteristic function 
rewrites 

^n 3 (p)-i(M p (^)-E[M p ( /t )])W 

= expf /" ^ G (en 3 (p)- 1 fi(B(x,pr)))X(p)dxF(dr] 

\jR d xR+ 

Let n(dz) = 4>(z)dz with G L 1 ^) n L a (R d ). Since as p -> 
0n 3 (p)-V(S(x,pr)) - A(p)- 1 / a c d r<V(x) 
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and X(p) — > +00 

lim A(p)* G (en 3 (p)- 1 KB(x,pr))) 

-a a c2\0\ a r ad \<t>(x)\<* (l - e{6cp(x)) tan (™) b) 



dx almost everywhere, and this latter function is integrable with respect to dxF{dr) 
since (f> G L a (R d ) and J* R+ r ad F(dr) < +00. Furthermore, with 4>* given in (|27l) . we 
derive the following bound: 

|A(p)* G (en3(A»)-V(S(x,pr)))| < \(p)Cn 3 (p)- a HB(x,pr))\ a 

< Cr ad \4>*{x)\ a . 

This upper bound is independent of p and integrable with respect to dxF(dr) since 
(ft* G L a (M ). The dominated convergence theorem yields: 

lim / ^ G (9n 3 (p)- 1 p(B(x,pr)))X(p)dxF(dr) 

= -a a c2\9\ a [ r ad F(dr) I \(f>(x)\ a (l - e(9<f>(x)) tan{ira/2)b)dx. 

JU+ JR d 

This proves Theorem 12.191 □ 
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